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摘 要 : 基于 Shuttleworth-Wallace Hu(SWH) 双 源 蒸 散 模型 对 青藏 高 愿 那曲、 纳木错 、 藏 东南 站 蒸 散 
发 进行 估算 ,在 结果 验证 良好 基础 上 ,对 青藏 高 原 蒸 散发 变化 特征 及 各 站 主要 影响 因素 进行 了 分 
析 。 结 果 表 明 :SWH 模 型 在 青藏 高 原 3 个 草 旬 站 适用 性 良好 ;年 茹 散发 介 于 388~732 mm 之 间 , 年 
内 分 布 呈 先 增 大 后 减 小 特征 ;3 站 燕 散 发 组 分 差异 较 大 ,那曲 站 和 纳木错 站 土壤 燕 发 对 燕 散 总 量 的 
贡献 分 别 为 53% 和 56%, 藏 东南 站 蒸 散发 则 几乎 全 部 由 植被 莱 腾 贡献 , 占 比 高 达 9596; 植 被 叶 面积 
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青藏 高 原 地 处 中 国 西南 部 ,是 世界 上 海拔 最 
高 .面积 最 大 、 地 形 最 为 复杂 的 高 原 , 被 称 为 “世界 
屋脊 “ 址 界 第 三 极 ”' 。 青 藏 高 原 63.5% 以 上 的 面 
积 被 高 山 荒漠 草原 ,高山 草原 和 草 多 覆盖 。 随 着 全 
球 变 暧 ,以 青藏 高 原 为 突出 代表 的 全 球 气候 和 环境 
都 发 生 了 显著 变化 ,高 原水 资源 呈现 持续 减少 趋 
势 ,水 循环 的 急剧 变化 深刻 影响 着 高 原 及 其 周边 地 
区 水 分 和 热量 平衡 一 。 

78 BUR. (Evapotranspiration ,ET) 是 地 表土 壤 蒸 
发 和 植被 蒸腾 的 总 和 。 蒸 散发 是 陆 表 - 大 气 间 进 行 
水 -能 量 交 换 的 主要 过 程 , 在 陆地 水 循环 、 碳 循环 以 
及 能 量 平衡 过 程 中 起 着 决定 性 作用 “ 。 在 水 循环 
过 程 中 ,陆地 上 大 约 67% 的 降水 会 通过 蔡 散 返回 大 
气 中 中 。 生 态 系 统 燕 散 发 主要 包括 土壤 蒸发 和 植被 
入 认识 水 循环 过 程 的 重要 途径 。 目 前 许多 研究 者 利 
JH Penman-Monteith 模型 .能 量 平 衡 模 型 ”和 遥感 
数据 产品 “5 等 估算 蒸 散发 , 均 取 得 了 较 理 想 的 结 
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果 。 利 用 传统 模型 估算 蒸 散 发 受到 气象 站 点 分 布 
不 均 的 影响 ,难以 对 大 范围 区 域 尺 度 的 蒸 散 发 进行 
(HE: ; 而 利用 GLEAM GLDAS, MODIS16 45: 3 RAE 
散 产 品 虽 然 解决 了 大 范围 区 域 斥 度 蒸 散发 的 估算 
问题 “” ,但 难以 拆 分 蒸 散 发 组 分 并 对 其 进行 定量 
分 析 o Shuttleworth 和 Wallace 在 1985 年 提出 了 
Shuttleworth-Wallace(S-W ) 双 源 蒸 散 模 型 2 ,实现 了 
对 薰 散 总 量 进行 估算 的 同时 又 对 其 组 分 进行 拆 分 。 
S-W 模 型 是 一 个 物理 机 制 明确 的 燕 散发 概念 模型 ,经 
阻抗 参数 化 以 后 ,具备 很 强 的 土壤 蒸发 -植被 蒸腾 模 
拟 能 力 ,被 广泛 应 用 于 莹 散发 估算 研究 7 s Kato 
等 和 Brisson 等 利用 实测 数据 对 S-W 模型 模拟 
蒸 散发 结果 进行 验证 ,分 析 了 实测 蒸 散 发 与 模拟 薰 
散发 之 间 的 差异 ,结果 均 表 明 S- 双 模型 在 稀 玻 冠 层 
蒸 散 发 估算 中 具有 较 高 的 准确 性 。 

Hu 47° 通过 引入 Ball-Berry 气孔 导 度 模型 和 
基于 光 能 利用 率 的 植被 总 初级 生产 力 (Gross prima- 
ry productivity, CPP) 模 型 ,实现 了 基于 气象 数据 和 
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GPP 对 冠 层 气 孔 导 度 的 估算 ,形成 了 新 的 SWH 双 源 
菩 散 模型 。Hu 等 ”将 SWH 模 型 应 用 在 中 国 长 白山 
站 (森林 生态 系统 ) 和 海北 站 (草地 生态 系统 ) ,取得 
了 较为 理想 的 模拟 效果 。 吴 戈 男 等 "在 全 国 51 个 
站 点 对 SWH 模 型 模拟 结果 与 各 站 多 年 蒸 散 观 测 数 
据 在 年 和 季节 尺度 上 进行 了 对 比 验 证 ,结果 表明 
SWH 模型 具有 良好 的 模拟 效果 。Jiang 等 所 利用 
SWH 模型 对 1981 一 2010 年 黄河 流域 蒸 散发 进行 了 
估算 并 对 其 组 分 进行 了 拆 分 。 上 述 研究 表明 ,SWH 
模型 能 够 较 好 地 模拟 站 点 和 区 域 尺度 陆 面 蒸 散 ,并 
且 能 够 对 蒸 散 发 组 分 进行 拆 分 ,这 为 深入 分 析 生 态 
系统 蒸 散 全 过 程 提 供 了 途径 。 

青藏 高 原 也 被 称 为 “亚洲 水 塔 ”, 是 亚洲 几 条 主 
要 河流 的 源头 。 青 藏 高 原 的 降水 和 蒸 散 发 决定 着 
下 游 超过 2x10" 人 的 淡水 供应 。 作 为 水 分 和 能 量 : 
衡 的 重要 环节 ,青藏 高 原 蒸 散 及 其 变化 特征 受到 广 
泛 关注 ,然而 如 何 细 化 蒸 散 发 过 程 仍 是 目前 研究 的 
重点 。 此 外 ,量化 青藏 高 原 蒸 散 发 . 拆 分 其 组 分 对 
水 循环 过 程 研究 及 认识 高 原 区 域 生 态 系统 .应 对 气 
候 变 化 国家 战略 的 制定 乃至 国家 安全 具有 重要 的 
战略 意义 。 本 研究 评估 了 SWH 模 型 在 青藏 高 原 地 
区 模拟 蒸 散 发 的 精度 ,并 基于 结果 分 析 蒸 散发 及 其 
组 分 变化 特征 ,探讨 影响 蒸 散 发 变化 的 主要 因素 ， 
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注 :该 图 基于 国家 测绘 地 理 信 息 局 标准 地 图 服务 网 站 下 载 的 
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为 后 续 青 藏 高 原 区 域 尺 度 上 燕 散发 估算 及 燕 散 发 
组 分 时 空格 局 与 变化 规律 等 研究 提供 科学 依据 。 


1 研究 区 概况 


本 文 所 选 青藏 高 原 3 个 高 寒 草 向 观测 站 位 置 如 
图 1。 那 曲 高 寒气 候 环境 观测 研究 站 海拔 4509 m, 
地 处 唐古拉 山脉 与 念 青 唐古拉 山脉 之 间 的 那曲 河 
流域 ,位 于 羌 塘 高 原 的 东部 和 怒江 流域 的 上 游 ,分 
布 有 季节 性 和 多 年 冻 土 。 那 曲 站 所 处 位 置 植被 覆 
盖 良 好 , 草 密集 ,土壤 主要 为 砂 质 粉 砂 壤 十 WA 
错 多 圈 层 综合 观测 研究 站 海拔 4730 m ,位 于 纳木错 
湖岸 ,其 后 是 念 青 唐 古 拉 山 。 纳 木 错 站 被 高 山 草 多 
覆盖 ,土壤 类 型 主要 为 砂 质 粉 砂 塘 土 ,砾石 含量 在 
地 面 以 下 30~40 cm 处 较 高 。 藏 东南 高 山 环境 综合 
观测 研究 站 位 于 靠近 森林 覆盖 的 青藏 高 原 东 南部 
的 山谷 中 ,海拔 3327 m ,植被 成 密 , 主要 是 高 寒 草 向 
和 温带 针 叶 林 。 


2 数据 与 方法 


2.1 数据 来 源 
本 文选 取 那 曲 .纳木错 、 藏 东南 站 2006—2010 
年 小 时 尺度 连续 观测 气象 环境 要 素数 据 ,数据 来 源 


À 


图 号 为 GS(2019)1822 号 的 标准 地 图 制作 , 底 图 边界 无 修改 。 


图 1 站 点 位 置 分 布 


Fig. 1 Location distribution of the stations 
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于 国家 青藏 高 原 科学 数据 中 心 (http:/data.tpdc.ac. 
cn)", fign f182& (Leaf area index ,LAT) 数 据 
来 源 于 中 国 科学 院 资源 环境 科学 与 数据 中 心 (https: 
//www.resde.cn ) $2 FE AY GLOBMAP 全 球 叶 面积 指数 
产品 (Version 3) ,该 产品 基于 AVHRR il MODIS 数 
据 定量 融合 反 演 得 到 ,空间 分 辩 率 为 8 km, 时 间 分 
辩 率 为 8 d。 由 于 LAI 数 据 为 8 4 时间 尺度 ,进一步 
将 模型 驱动 数据 进行 处 理 , 最 终 模型 所 有 输入 和 输 
出 数据 均 为 8 d 平 均值 。 

本 文采 用 中 国 陆地 ET 产品 数据 (GPR) 和 青藏 
高 原 陆 表 月 蒸 散发 数据 (ET.) 对 SWH 模 型 模拟 蒸 散 
发 进行 验证 评价 ,通过 对 模型 模拟 蒸 散发 结果 和 蒸 
散发 产品 数据 进行 线性 回归 分 析 ,采用 决定 系数 
(R) .平均 绝对 误差 (MAE ) L907 TRO 23 (RMSE) LA 
及 平均 相对 误差 (MRE)4 种 评价 指标 来 检验 模型 模 
TAZSBUR BE. GPR 数据 产品 包含 2000 一 2018 年 
中 国 陆 表 蔡 散发 ,空间 分 辨 率 为 1 km, 时 间 分 辩 率 
为 10 d;ET, 数 据 产 品 包含 2000—2018 年 青藏 高 原 
月 平均 陆 表 蒸 散 发 ,空间 分 辨 率 为 0.01*x0.01” ,时 
间 分 辩 率 为 月 尺度 。 本 文 所 用 验证 数据 均 来 源 于 
国家 青藏 高 原 科 学 数据 中 心 (http://data.tpdc.ac.cn)。 
2.2 SWH Ate! 

Hu 等 “通过 引入 土壤 表面 阻抗 方程 、Ball-Ber- 
ry 气孔 导 度 模型 和 基于 光 能 利用 率 的 GPP 模 型 ,在 
S-W 模 型 的 基础 上 形成 了 新 的 SWH 模 型 。 经 阻抗 
参数 化 、 冠 层 气 孔 导 度 估算 等 改进 ,SWH 模 型 具备 
了 很 强 的 土壤 蒸发 和 植被 蒸腾 模拟 能 力 ,能 够 针对 
不 同 生态 系统 进行 关键 参数 取 值 ,从 而 实现 在 不 同 
生态 系统 中 的 良好 模拟 性 能 ,可 广泛 应 用 于 各 类 生 


态 系统 蒸 散 估算 ”” 。 
模型 核心 公式 如 下 : 
AET=AT+AE =C PM, +C PM, (1) 
pc, VPD- Ar, R, 
P 
PM, = um ac (2 ) 
5. 
A i 5 A laa + 一 
"Tm pc, VPD - Bratt ER.) 
PM, = Tu Tas (3 ) 
I. 
a à 5 2 Tu + -| 
C = (4) 
PPa 


C=— 1 (5) 
: P.P, 


pp. *p.) 

式 中 :ET 为 模拟 蒸 散 发 (mm) ;7 为 植被 蒸腾 (mm ) ; 
E 23 ABER (mm); A 为 液态 水 蒸发 为 水 汽 的 相 变 
HAMI- kg) ; PM, WER E E (mm); PM, HE 
IRR Az 5t (mm) ; C, 为 植被 蒸腾 比例 系数 ; C. 为 土 
SEAR A AA BL; A 为 饱和 水 汽 压 随 温度 变化 曲线 
HIRR (KPa K); R 为 冠 层 项 部 获得 能 量 (MJ.m”); 
p 为 空气 密度 (kg"m ); e, 为 空气 定 压 比 热 容 ; R, 
为 土壤 表层 吸收 能 量 (MJ.m2); VPD 为 饱和 水 汽 压 
差 (kPa); y 为 湿度 计 常 数 (kPa'K2); p, 、p.、p. 为 
中 间 变 量 ; ro ru ru ru 和 六 分 别 为 地 表 与 冠 层 
高 度 间 的 空气 动力 学 阻抗 . 冠 层 高 度 与 参考 高 度 间 
的 空气 动力 学 阻抗 \ 冠 层 边界 层 阻抗 . 冠 层 气孔 阻 
抗 和 土壤 表面 阻抗 (sm )。 

SWH 模 型 驱动 数据 包括 气象 环境 要 素数 据 与 
遥感 数据 ,气象 环境 要 素数 据 包括 平均 气温 ( 思 ) 、 相 
对 湿度 (RH) 饱和 水 汽 压 差 (VPD) KE (Wind) 24 
辐射 (及 ) .土壤 含水 量 (SW) .土壤 热 通 量 (C) .光合 
有 效 辐 射 (PAR) ,遥感 数据 为 LAI 数 据 。 通 过 驱动 
数据 求 出 上 > 、r rr sr, 5 种 阻抗 ,并 进一步 计 
算 土 壤 蒸发 和 植被 蒸腾 ,从 而 求 得 蒸 散 发 。SWH 模 
型 详细 介绍 参见 Hu 等 ”” ,本 文 模型 模拟 过 程 中 重 
HPUH Hu E ba bs ai 2o En d BURAYI 
为 3.4.508 ,10.3 ,0.017 ,0.0012 .188。 


3 结果 与 分 析 


3.1 SWH 模 型 模拟 评估 

模型 主要 输入 数据 2006 一 2010 年 8d 变 化 情况 
如 图 2。 那 曲 纳木错、 藏 东南 站 各 要 素 变 化 大 致 相 
同 ,整体 均 呈 现 随时 间 变 化 的 年 内 和 年 际 变化 ,但 
个 别 要 素 差 异 较 大 。 纳 木 错 站 VPD 整体 高 于 那曲 
站 和 藏 东南 站 ;那曲 站 SW 明显 高 于 纳木错 站 和 藏 
东南 站 ,其 中 藏 东南 站 SW 最 小 ;R, 和 G6 变化 中 , 藏 东 
南 站 波动 最 为 明显 ,其 余 2 站 变化 幅度 较 小 ;3 站 
LAI 变 化 差异 较 大 , 藏 东南 站 日 均 LAI 达 30, 而 那曲 
站 和 纳木错 站 仪 为 2.4 和 0.9。 

SWH 模 型 模拟 验证 结果 如 图 3。 青 藏 高 原 3 站 
模拟 ET 月 值 与 GPR 数据 、ET, 数 据 拟 合 结果 均 较 
好 ,那曲 .纳木错 、 藏 东南 站 与 GPR 拟 合 尼 分 别 为 
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(c) RH 
一 那曲 站 ---- 纳木错 站 
一 - 藏 东南 站 


Pa， d 
i 


RH/% 
2 
= 


e «ww d (f 
vs $9. S $ 
日 期 (年 -月 -日 ) 日 期 (年 -月 -日 ) 日 期 (年 -月 -日 ) 
(d) VPD (e) SW (f) Ry 
一 那曲 站 ---- 纳木错 站 一 那曲 站 ---- 纳木错 站 那曲 站 ---- 纳木错 站 
"EET 35 。 一 藏 东南 站 30, ~~ 藏 东南 站 
st * 
E. x t 
a > ; 
& [2] 
3 
L3 
s Xe us Ss 
日 期 (年 -月 -日 ) 日 期 (年 -月 -日 ) 
(h) Wind (i) LAI 
一 那曲 站 ---- 纳木错 站 — 那曲 站 -…- 纳木错 站 
— 藏 东南 站 


G/MJ-m?-d" 
Wind/m-s* 


日 期 (年 -月 -日 ) 


日 期 (年 -月 -日 ) 


LAI 


日 期 (年 -月 -日 ) 


E: PAR 为 光合 有 效 辐射 ;了 ,为 了 


F 均 气温 ;RH 为 相对 湿度 ;VPD 为 饱和 水 汽 压 差 ;SW 为 土壤 含水 量 ;RR 为 净 辐 射 ;6 为 土壤 热 通 量 ; 
Wind 为 风速 ;LAI 为 植被 叶 面 积 指数 。 下 同 。 
图 2 3 站 SWH 模 型 输入 数据 日 变化 


Fig. 2 Diurnal variations of SWH model input data at three stations 


0.89 , 0.69, 0.69, RMSE 4} 5!) X 9.55 mm + d ', 13.10 
mm *d' £l 14.33 mm- d! , MAE 47-4) 4 7.75 mm- d^' , 
15.68 mm- d” ,16.94 mm- d” , MRE 分别 为 49% ,12596 
fll 5196 ; 4G ET, WE R^ 4) 31] 2g 0.88, 0.68 和 0.71, 
RMSE 分 别 为 10.10 mm* d^ ,13.25 mm- d^' Fil 13.75 
mm*d'', MAE 7 il Jj 9.18 mm* d^, 12.09 mm- d! , 
15.90 mm d^, MRE 分 别 为 29% , 53% FI 50%. V. 
SWH 模 型 模拟 ET 与 验证 数据 月 变化 情况 可 看 出 ， 
SWH 模 型 模拟 结果 与 验证 数据 变化 趋势 基本 一 致 
(图 4) ,相对 2 套 验 证 产品 而 言 ,SWH 模 型 模拟 ET 


虽然 存在 略微 低估 或 高 佑 的 现象 ,但 整体 变化 较为 
平滑 稳定 。GPR 产品 和 ET, 产 品 与 模型 模拟 结果 对 
比 精度 不 完全 一 致 ,这 可 能 与 2 套 产品 的 分 辩 率 、 观 
测 手 段 和 数据 融合 方法 等 存在 差异 有 关 。 通 过 验 
证 分 析 , 认 为 SWH 模 型 在 青藏 高 原 站 点 尺度 具有 较 
好 的 适用 性 ,可 进一步 用 于 青藏 高 原 3 个 草 甸 站 蒸 
散发 模拟 人 研究。 
3.2 不 同时 间 尺 度 蒸 散发 变化 

2006—2010 4E. 3 x H ET ZE fL Al 5,3 xh H z& 
散发 整体 呈 先 增加 后 减少 的 趋势 。2006 一 2010 年 
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注 :ET 为 蒸 散发 ;GPR 为 中 国 陆地 ET 产品 数据 ;ET, 为 青藏 高 原 陆 表 月 蒸 散 发 数据 。 下 同 。 


图 3 模拟 ET 与 GPR ET 对 比 


Fig. 3 Comparison of modeled ET with GPR and ET, 


(a) 那曲 站 
120 p  -— 模拟 ET —— GPR 
— ET, 


月 ET/mm 


日 期 (年 -月 ) 


(b) 纳木错 站 
120 -  —— 模拟 ET 一 - GPR 


日 期 (年 -月 ) 


图 4 模拟 ET 与 GPR .ET 月 变化 


(c) 藏 东南 站 
120 一 一 模拟 ET —- GPR 
—— ET, 


日 期 (年 -月 ) 


Fig. 4 Monthly variations in modeled ET, GPR, and ET, 


间 ,那曲 纳木错、 藏 东南 站 日 ET 分 别 介 于 0.0~3.3 
mm ,0.1~3.3 mm,0.3-4.9 mm,3 站 日 ET 均值 分 别 为 
1.2 mm,1.4 mm,1.9 mm。 总 的 来 看 ,3 站 日 ET 变化 
趋势 较为 相似 ,但 那曲 站 日 ET 变化 相对 平缓 , 较 纳 
木 错 站 和 藏 东南 站 波动 小 .稳定 性 高 , 藏 东南 站 变 
化 最 为 剧烈 ,纳木错 站 变化 波动 最 大 。 

3 站 蒸 散发 年 内 月 变化 如 图 6, 那 曲 ` 纳 木 错 i 


东南 站 月 ET 均 呈 先 增加 后 减少 的 单 峰 型 变化 趋势 ， 
但 各 站 变化 有 所 差异 。1 一 3 月 各 站 ET 均 较 小 且 较 
为 稳定 ,4 一 7 月 3 站 ET 均 有 不 同 程度 的 增加 ,其 中 ， 
7 月 那曲 站 (85 mm) 、 藏 东南 站 (97 mm) H ET isi 
大 ,纳木错 站 ET 最 大 值 出 现在 8 月 ,为 77 mm。 总 
体 上 , 藏 东南 站 月 ET 大 于 那曲 站 和 纳木错 站 。ET 
年 内 月 变化 与 植被 生长 季 (3 一 11 月 ) 基 本 吻合 ,1 一 
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静 等 ， 基 于 SWH 模 型 的 青藏 高 


5 一 一 那曲 站 


原 高 寒 草 多 蒸 散 发 及 其 组 分 变化 分 析 
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- 纳木错 站 


oc 藏 东南 站 : 


o ^ © ^ © ^ 
~ p^ ~ p^ ~ % 
SÉ QU SU GF Ux 
S S S g SS 
ty ty ty ty ty ty 
日 期 (年 -月 -日 ) 


5 3 站 日 ET 变化 
Fig. 5 Variations in daily ET of three stations 
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6 3 站 月 ET 变化 
Fig.6 Variations in monthly ET of three stations 


2 月 ET 变化 不 大 ,3 月 起 ET 增长 速度 加 快 ,特别 是 4 
月 之 后 3 站 ET 均 呈 现 明 显 的 增长 趋势 ,直到 达到 最 
大 后 减 小 ,11、12 月 ET 与 1.2 月 相近 。 另 外 ,高 原 季 
风 的 变化 会 给 区 域 气象 带 来 较 大 影响 ,高 原 季 风 一 般 
6 月 来 临 ,7 一 9 月 持续 加 强 , 到 10 月 逐渐 减退 ” ,季风 
会 使 降水 增加 , 净 辐 射 增强 ,从 而 使 得 高 原 蒸 散 发 增 
加 ,这 也 在 一 定 程度 上 影响 了 燕 散 发 月 变化 。 

图 7 为 3 站 莹 散发 年 及 季节 变化 ,3 站 蒸 散 发 年 
际 变化 不 大 ,整体 略微 上 升 或 不 变 。 那 曲 、 纳 木 错 、 
藏 东南 站 年 ET 分 别 介 于 388~479 mm ,450-658 mm 
和 674~732 mm 之 间 ,3 站 年 ET 均 在 2009 年 达到 最 


大 ,分 别 为 479 mm、658 mm 和 732 mm, ZN EE 
上 ,那曲 站 ET 四 季 变 化 波动 不 大 ,其 中 冬季 几乎 没 
有 变化 ;纳木错 站 IT 四 季 变 化 均 不 稳定 ,其 中 冬季 
波动 最 明显 .最 剧烈 ; 藏 东南 站 ET 春 、 夏 季 变 化 均 有 
缓慢 上 升 的 趋势 , 秋 冬季 出 现 了 较为 明显 的 波动 。 
3.3 蒸 散 发 组 分 分 离 

通过 SWH 模 型 拆 分 得 到 那曲 、 纳 木 错 、 藏 东南 
站 2006 一 2010 年 土壤 蒸发 有) 和 植被 蒸腾 (7)。 图 
8a~c 可 知 ,那曲 站 日 五 变化 范围 为 0.0~1.1mm, 日 了 
介 于 0.0~2.9 mm 之 间 ,植被 非 生长 季 内 日 已 和 日 了 
变化 趋势 较为 吻合 ,生长 季 内 两 者 变化 呈 相 反 趋 
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图 7 3 站 年 和 季 尺 度 ET 变 化 


Fig.7 Variations in ET of three stations at interannual and seasonal scales 


势 , 即 7 不 断 增 加 ,E 则 呈 减 少 趋势 ;那曲 站 日 土壤 
蒸发 占 蒸 散 总 量 的 比值 CBET) 介 于 0.07~0.88 ,均值 
约 为 0.53。 纳 木 错 站 日 变化 波动 大 ,最 大 值 约 为 
1.3 mm, 最 小 值 接近 于 0; 纳 木 错 站 日 7 最 大 值 约 为 
2.5 mm; 纳木错 站 日 B/ET 波 动 较为 剧烈 ,最 大 值 约 
为 0.89, 最 小 值 接近 于 0, 日 B/ET 均 值 为 0.56。 藏 东 
南 站 日 均值 约 为 0.1 mm, A TURAE TRAE, 
且 变 化 趋势 与 日 ET 相似 ,日 ET 均值 为 0.05。3 站 
蔡 散 发 组 分 年 内 月 变化 如 图 8d~f。 那 曲 站 月 变化 
整体 较为 稳定 ,1 一 5 月 呈 小 幅 上 升 趋势 ,年 内 月 最 
大 值 出 现在 5 月 ,为 22 mm; 月 7 呈 单 峰 型 变化 ,1 一 5 
月 11 一 12 月 较为 稳定 ,7 月 (76 mm) 达 到 最 大 ;月 尺 
度 BEVET 呈 先 减少 后 增加 的 趋势 ,其 中 ,1 一 4 月 较为 
稳定 ,5 一 8 月 不 断 减 小 ,而 后 呈 增 长 趋势 ,最 小 值 出 
现在 8 月 ,为 0.09, 最 大 值 出 现 3 月 (0.82)。 纳 木 错 
站 月 7 在 1 一 4 月 .10 一 12 月 较为 稳定 ,5 一 10 月 呈 先 
增加 后 减 小 趋势 ,最 大 值 出 现在 8 月 ,为 49 mm; H E 


在 1 一 3 月 呈 上 升 趋势 ,而 后 基本 在 30 mm 左右 小 
幅 波 动 ;纳木错 站 E/ET 月 变化 呈现 2 个 峰值 ,分 别 
在 3 月 (0.84) 和 11 月 (0.81) ,最 小 值 出 现在 8 月 
(0.37)。 藏 东南 站 月 E 小 旦 稳定 ,月 均 E 为 2.4 mm; 
月 7 则 呈 先 增加 后 减 小 趋势 ,1 一 4 月 缓慢 上 升 , 变 化 
较 稳 定 ,4 一 7 月 上 升 趋势 明显 ,最 大 值 出 现在 7 月 ， 
为 95 mm,7 一 12 月 则 不 断 减少 ; 藏 东南 站 月 E/ET 整 
体 较 小 ,波动 范围 为 0.02~0.11, 月 E/ET 均 值 为 
0.05。 藏 东南 站 几乎 所 有 蒸 散 发 均 由 植被 蒸腾 页 
献 ,由 于 该 站 草 多 发 育 良 好 , 且 混 有 温带 针 叶 林 , 植 
被 生长 茂 感 ,生态 系统 蒸 散 主要 来 源 于 植被 蒸腾 。 
3.4 蒸 散发 影响 因素 分 析 

薰 散 发 主要 受气 温 .降水 .风速 日照 时 数 等 气 
象 要 素 和 植被 生长 状况 等 诸多 因素 影响 , 而 不 同 
素 之 间 也 相互 影响 , 薰 散发 的 变化 成 因 十 分 复杂 。 
通过 对 3 个 站 点 蒸 散 发 与 主要 驱动 因子 作 相 关 分 
Wr ,选取 与 燕 散发 决定 系数 ( 尼 ) 大 于 0.7 的 驱动 因子 
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(b) 纳木错 站 日 变化 
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注 :7 为 植被 蒸腾 ;五 为 土壤 蒸发 ;BET H ERRE UY LLL. 
图 8 日 和 月 尺度 ET 及 其 组 分 变化 


Fig.8 Variations in daily and monthly ET and its components 


进一步 辨析 影响 蒸 散 发 的 关键 要 素 。 采 用 分 层 回 
归 模 型 检验 主要 影响 因素 对 3 站 莹 散发 影响 程度 及 
重要 性 ( 表 1)。3 站 仅 考 虑 单一 要 素 影响 时 ,结果 并 
不 理想 , 当 综合 纳入 主要 影响 因素 后 , 尼 均 超过 
0.8。 从 标准 化 系数 可 知 ,LAI 对 3 站 蒸 散 发 影响 均 
较 其 他 因素 大 ,VPD 在 藏 东南 站 对 蒸 散发 的 影响 仅 
次 于 LAI, 以 上 结果 均 通 过 0.001 显著 性 水 平 检验 。 
LAI 与 植被 生长 状况 密切 相关 ,其 通过 影响 植 
被 蒸腾 影响 蒸 散发 过 程 。 以 往 研 究 中 ,学 者 们 就 
VPD 对 莹 散发 的 影响 持 不 同 态度 ,Dai 8$ 34A DI eg 
VPD 值 就 像 空气 中 的 “水 和 泵 ”, 从 地 表土 壤 和 植被 中 
吸收 水 分 ,增加 了 裸露 土壤 的 水 分 流失 。 植 被 长 期 
处 于 水 分 荐 乏 条件 时 ,VPD 增加 可 能 会 导致 气孔 导 


度 降低 ,抑制 植被 生长 ,降低 植被 蒸腾 作用 。 而 
Grossiord 等 3 使 用 9 个 气孔 导 度 模型 研究 VPD 如 何 
影响 蒸腾 作用 ,发 现 VPD 增 大 导致 气孔 导 度 降低 ， 
但 在 给 定 的 VPD 阔 值 下 蒸腾 作用 增强 了 。 张 亚 春 
等 % 通过 分 析 青 藏 高 原 不 同 下 垫 面 的 5 个 观测 站 点 
蒸 散发 影响 因素 发 现 , 实 际 日 莱 散 发 与 气温 .地温 、 
净 辐 射 和 土壤 含水 量 的 相关 性 较 大 , 且 均 呈 显 著 正 
相关 。Ma 等 认为 降水 是 控制 青藏 高 原 大 部 分 地 
区 ET 变化 的 主要 驱动 因素 ,Wang 等 “认为 LAL 和 
降水 是 ET 的 主导 变量 ,而 土壤 蒸发 则 由 VPD .平均 
气温 \ 净 辐射 和 降水 控制 。 本 文 结论 与 上 述 研 究 存 
在 相似 和 差异 之 处 ,可 能 与 研究 范围 驱动 数据 Ex 
型 方法 等 不 同 有 关 。 
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表 1 蒸 散 发 与 各 要 素 的 分 层 回 归 及 未 标准 化 系数 (标准 化 系数 ) 
Tab. 1 Stratified regression of ET with each element 
站 点 要 素 模型 1 模型 2 模型 3 模型 4 模型 5 模型 6 
那曲 站 PAR 0.085™(0.726) 0.002(0.019) -0.000(0.000) -0.017™*(-0.146) -0.011(-0.096) — -0.005(-0.042) 
T. i 0.106"(0.920) | 0.1057(0.909) 0.075"(0.655) | 0.0737(0.636) 0.005 (0.044) 
VPD - - 0.229(0.035) 0.648 (0.097 ) 0.9057(0.136) —— 1.7197 (0.259) 
R, 一 - - 0.0827(0.348) | 0.0027(0.390) — 0.085" (0.359) 
G - - - - -0.272'(-0.115) 0.068 (0.029) 
LAI - - - - - 0.190“(0.485 ) 
Constant -1.369 1.158 1.153 0.965 0.634 -0.237 
N 230 230 230 230 230 230 
R 0.528 0.875 0.875 0.888 0.891 0.952 
AR 0.528 0.347 0.000 0.013 0.002 0.061 
纳木错 站 T 0.0927 (0.835) 0.066™(0.599) 0.035"(0.320) - - - 
R, - 0.054"(0.328) | 0.0477(0.289) - - - 
LAI - - 0.590°"(0.405) - - - 
Constant 1.495 1.123 0.597 - 一 一 
N 230 230 230 - - - 
R 0.697 0.748 0.818 - - - 
AR 0.697 0.052 0.069 - - - 
藏 东南 站 PAR 0.1177(0.787) 0.074™(0.497) 0.048 (0.320) 0.028 (0.190  - - 
T. - 0.066 (0.440) 0.0547(0.359) -0.008(-0.056) - - 
VPD - - 4.196""(0.326) 5.092"(0.396)  - - 
LAI - - - 0.043°"(0.496)  - - 
Constant -1.177 —0.388 —0.589 -1.275 一 m 
N 138 138 138 138 - - 
R 0.619 0.729 0.779 0.807 - - 
AR’ 0.619 0.110 0.050 0.029 一 


注 :* 4) FI EAN P«0.05 .P<0.01、P<0.001。PAR 为 光合 有 效 辐 射 ;7 为 平均 气温 ;VPD 为 饱和 水 汽 压 差 ;有 ,为 净 辐 射 ; 6 为 土壤 热 通 量 ; 


LAI 为 植被 叶 面积 指数 ;Constant 为 常数 ;为 样本 数 ; 尼 为 决定 系数 ;AP 为 决定 系数 变化 量 。 


4 讨论 


以 往 研 究 认 为 植被 蒸腾 是 全 球 蒸 散发 的 主体 
部 分 ,但 Kool 等 号 证 明 土 壤 蒸发 是 干旱 和 半 干 旱地 
区 蒸 散 发 重要 的 组 成 部 分 ,Ma 5e 538 et PML 模型 
计算 得 出 青藏 高 原 蒸 散发 64% 来 源 于 土壤 蒸发 ， 
Wang 等 中 同样 利用 PML 模 型 估算 了 青藏 高 原 实 际 
蒸 散 发 并 对 蔡 散 发 组 分 进行 了 量化 ,发 现 除 位 于 青 
藏 高 原 东 南部 具有 高 树冠 覆盖 的 半 湿 润 地 区 外 , 土 
坏 营 发 均 为 蒸 散发 主要 组 成 ,上 述 研究 结果 与 本 文 
蒸 散 发 组 分 拆 分 结果 相似 。SWH 模 型 由 于 受到 通 
感 数 据 时 间 尺 度 的 影响 ,输出 结果 为 8d 均 值 数据 ， 
不 能 细 化 蒸 散 发 日 变化 ,对 本 文 燕 散发 日 尺度 变化 
分 析 带 来 了 一 定 程度 的 误差 ,Zhao HK METRIC 


模型 与 SWH 模 型 结合 ,弥补 了 SWH 模 型 在 日 尺度 
上 模拟 蒸 散发 的 精度 缺陷 。 本 文 3 个 站 点 蒸 散发 及 
其 组 分 变化 有 较 大 差异 ,3 个 观测 站 虽 均 为 高 寒 草 
多 下 热 面 ,但 各 站 环境 、 气 候 条 件 等 有 较 大 差异 。 
受到 河流 湖泊 的 影响 ,那曲 站 和 纳木错 站 蒸 散发 较 
受到 高 寒 草 甸 和 温带 针 叶 林 同时 覆盖 影响 的 藏 东 
南 站 小 ,而 本 文 在 模拟 3 个 站 点 蒸 散发 时 重要 参数 
取 值 为 草地 生态 系统 经 验 参 数 ,也 使 得 模拟 结果 可 
能 与 蒸 散发 实际 变化 有 所 差异 。 后 续 要 将 SW 瑞 模 
型 应 用 于 整个 青藏 高 原 地 区 ,如 何 解决 区 域内 环 
境 .气候 等 对 模型 结果 造成 的 影响 仍 需 深 入 思 

SWH 模 型 只 对 土壤 蒸发 和 植被 蒸腾 进行 了 拆 分 ,未 
对 植被 冠 层 截留 蒸发 进行 估算 ,后 续 也 需 将 冠 层 蒸 
发 纳 和 人 蒸 散发 拆 分 人 研究 ,以 期 更 为 准确 地 描述 青藏 
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Variations of evapotranspiration and its components in alpine meadow 
on the Tibetan Plateau based on SWH model 


MEI Jing, SUN Meiping", LI Lin' 
(1. College of Geography and Environmental Science, Northwest Normal University, Lanzhou 730070, Gansu, China; 


2. Northwest Institute of Eco-Environment and Resources, CAS, Lanzhou 730000, Gansu, China) 


Abstract: Based on the SWH dual source model, the evapotranspiration of the Nagqu Station of Plateau Climate 
and Environment, the Nam Co Monitoring and Research Station for Multisphere Interactions, and the Southeast 
Tibet Observation and Research Station for the Alpine Environment was estimated. Based on the validation of the 
results, the characteristics of evapotranspiration changes on the Qinghai-Tibet Plateau, China and the main influencing 
factors of each station were analyzed. The main conclusions are as follows: the SWH model is well applicable to 
the three stations on the Tibetan Plateau. The annual evapotranspiration ranges from 388 mm to 732 mm, with an 
increase in interannual distribution, followed by a decrease. The differences in evapotranspiration components 
between the three stations are significant, the contribution of evaporation to evapotranspiration at the Nagqu and 
Nam Co stations is 53% and 56%, respectively, and the evapotranspiration at the Southeast Tibet Station is almost 
entirely contributed by transpiration, accounting for 95%. The leaf area index is the most important factor 
affecting evapotranspiration at the three stations. The vapor pressure deficit also has a significant effect on 
evapotranspiration at the Southeast Tibet Station. The findings of this study can provide a scientific basis for 
studying the temporal and spatial patterns of evapotranspiration and its components and the water cycle process 
over the Tibetan Plateau. 


Key words: evapotranspiration; components of evapotranspiration; SWH model; Qinghai-Tibet Plateau 


